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Abstract 


in  v'S*'f/jtrhed 

-^.'hts  reportAdeoumen‘t,.s  an  Investigation  --e^the  feasi¬ 
bility  of  usiqg  an  analytical  approach  and  the  vortex 
L(V±H>J  r+D  cr 

lattice  methodj^f of^ie  ^evaluatflMi  o-f  the  aerodynamic 

interference  effects  present  during  aerial  refueling. 

-  ?  A  KC-10  tanker  and  a  B-52  receiver  were  selected  for — - 

stud^^oweverV  the  method  oan  applies  to  any  tanker^0 

o 

receiver  combination.  Slxe^^major  assumptions  inolude: 
linear  potential  flow;  tanker  and  receiver  are  represented 
by  wing  planforms,’  fuselage  effects  areymall;  and  the 
rolling-up  process  of  the  tanke^  wing  tip  vortices  are 
not  considered.  The  analytical  approach  uses  a  lifting 
line  followed  by  a  semi-infinite  vortex  sheet  to  represent 
the  tanker.  Three  linear  lift  distributions  for  the 
tanker's  wing  were  used  and  equations  were  derived  for 


the  induced  downwash  at  any  point  in  apace.  Points  were 
aelectedr  on  the  receivers  wing\to  indicate  the  tanker's 
effect  on  the  receiver^  flowfield.  Results  of  the  analy- 

gj 

tical  equations  were  compared  with  analytical  predictions 
ajfrd?wer*  found  to  over^redict  the  induced  downwash  by 
25-35%*  ^iie  vortex  latex  lattiee-  eethed  used 

to  study  the  effects  of:  rectangular  vs  swept  wings  on 
induced  downwash}  tanker  tailplane  on  the  receiver;  change 
of  angle  of  attack  of  one  aircraft  due  to  the  presence  of 
the  other}  and  the  presence  of  the  tanker  on  the  receiver's 
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pitohing  moment.  *  Again*  the  results  for  induced  dovnvash 

were  oompared^ith  the  Douglas  predictions  to  evaluate  the 

aoouraoy.^VLM  results  varied  from  the  Douglas  prediction 

for  induced  dovnvash  by  only  5 %  and  were  even  more  aoourate 

. /  I 

at  lover  tanker  angles  of  attaok.*  Other  findings  ve.et 
swept  wings  were  necessary  for  good  simulation*  the  tanker 
tailplane  has  only  a  slight  effect  on  the  receiver*  and 
the  tanker  (in  this  combination  of  aircraft)  affects  the 
receiver  to  a  greater  extent  than  the  reoeiver  affects 
the  tanker.  The  results  indioate  that  a  simple  vortex 
lattice  method  which  ignores  fuselage  effeots*  can  be 
used  to  determine  many  of  the  interference  effects  exper- 
ienced  during  aerial  refueling. 


A  Stud*,-  of  the  Aerodynamic  Interference 
Effects  During  Aerial  Refueling 

I.  Introduction  . 

Background 

In  April  1923,  the  first  air  refueling  took  place 
between  two  Army  Air  Corps  DH-4B  aircraft.  Since  then  a 
variety  of  aircraft  has  been  used  as  both  tankers  and 
receivers.  In  the  past,  refueling  training  was  accom¬ 
plished  in  flight.  However,  with  the  rising  cost  of  both 
fuel  and  airframe  time,  the  Air  Force  became  interested 
in  decreasing  the  actual  flying  time  while  maintaining 
the  same  flying  proficiency.  This  was  especially  true 
with  heavy  (B-52,  C-141 »  C-5A)  type  aircraft.  The  most 
logical  way  to  accomplish  this  task  was  to  develop  a 
simulator  that  can  be  used  for  air  refueling  training  and 
pilot  proficiency  requirements. 

The  Simulator  Systems  Program  Office  of  Aeronautical 
Systems  Division  (ASD)  needs  a  method  to  provide  flow 
field  data  for  the  ^C-10  tanker  to  use  in  receiver 
simulator  design.  This  information  is  necessary  to  ensure 
that  the  receiver  simulator  "flies"  like  the  rtial  aircraft. 

During  aerial  refueling,  the  normal  position  of  the 
receiver  aircraft  is  below  and  behind  the  tanker.  In  the 
refueling  position,  the  presence  of  the  tanker  causes 
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significant  changes  in  the  aerodynamic  behavior  ox'  the 
receiver.  Therefore,  for  effective  refueling  simulation, 
it  is  of  prime  importance  to  determine  how  the  tanker 
aerodynamically  affects  the  receiver. 

Purpose 

The  objective  of  this  study  was  to  determine  if  an 
analytical  or  a  simple  numerical  approach  can  be  used  for 
the  evaluation  of  air  refueling  aerodynamics.  To  accomplish 
this,  it  was  necessary  to  identify  the  important  aerody¬ 
namic  factors  and  also  determine  the  complexity  of  the 
model  required  to  give  acceptable  results. 

The  following  effects  were  considered: 

1 .  Tanker  span-wise  lift  distribution  on 
the  flow  field. 

2.  Sweepback  of  the  tanker  and  the  receiver 
on  the  flow  field. 

3.  Tanker  tailplane  on  the  receiver. 

4*  Change  of  the  tanker  angle  of  attack 
due  to  the  receiver. 

5.  Change  of  the  receiver  angle  of  attack 
due  to  the  tanker. 

6.  Change  of  the  receiver's  pitching  moment 
due  to  the  tanker. 

The  fuselage  effects  of  both  tanker  and  receiver 
aircraft  were  not  considered  in  the  analysis.  This 
greatly  simplified  the  problem  and  comparison  of  the 
computed  values  with  the  available  Douglas  data  (5)» 
Indicated  that  it  was  a  reasonable  assumption.  A  dis¬ 
cussion  of  the  method  used  by  Douglas,  and  the  effects 
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they  considered,  is  found  in  Chapter  five. 

Approach 

The  study  of  the  aerodynamic  effects  of  one  aircraft 
on  another  is  not  a  new  problem.  In  fact,  aerial  refuel¬ 
ing  aerodynamics  is  similar  to  biplane  aerodynamics.  In 
biplane  theory,  the  vertical  distance  between  the  two  wings 
is  called  gap  and  the  horizontal  distance  between  the 
quarter  chords  of  the  two  wings  is  called  stagger.  The 
most  basic  method  to  model  the  tanker  and  receiver  is  to 
replace  each  aircraft  with  a  rectangular  wing.  The  re¬ 
fueling  situation  is  now  reduced  to  a  biplane'  problem 
with  an  exaggerated  stagger  and  gap. 

Glauert  (8),  Von  Mises  (18),  Munk  (13) ,  Millikan,  (12) 
and  others  developed  biplane  theory.  They  found  that,  for 
an  elliptical  lift  distribution  on  the  lead  wing,  a  numer¬ 
ical  solution  is  required  to  calculate  the  induced  normal 
velocity  at  any  point  in  the  flow  (10).  To  keep  the 
problem  analytical  and  identify  key  parameters,  three 
different  linear  distributions  for  the  tanker  wing  were 
selected.  They  are  the  tent  distribution  (Fig  l),  the 
modified  linear  distribution  (Fig  2),  and  a  single 
horseshoe  vortex  (Fig  3). 

To  extend  the  analytical  approach  and  investigate 
the  effects  of  sveepback  and  the  tanker  tailplane,  a 
vortex  lattice  method  (VLM)  was  used.  In  this  method, 
aircraft  wings  and  tailplanes  are  modeled  by  flat  plate 
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planforms  for  a  variety  of  configurations  and  positions. 

A  KC-10  tanker  and  a  B-52  receiver  were  used  in  the 
evaluation ,  however,  the  method  can  be  applied  to  any 
tanker  and  any  receiver. 

Result  Summary 

The  results  of  the  analytical  approach  and  the  vortex 
lattice  method  were  compared  wdth  Douglas  Aircraft  Company 
(Douglas)  data  (5).  It  was  found  that  the  analytical 
methods  were  in  qualitative  agreement  for  tanker  induced 
velocity,  while  the  vortex  lattice  results  showed  excel¬ 
lent  agreement  with  the  Douglas  data,  especially  at  lower 
angles  of  attack,  where  linear  theory  is  expected  to  give 
best  results. 


j 
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II 


•  Air  Refueling  -  A  Receiver  Pilot1  s  Perspective 

To  have  an  accurate  simulation,  it  is  necessary  to 
have  a  thorough  knowledge  of  the  process  that  is  simulated. 
The  same  is  true  for  the  refueling  problem.  This  chapter 
qualitatively  describes  the  air  refueling  process  from  a 
pilot's  point  of  view.  A  typical  refueling  of  a  KC-135 
tanker  and  a  B-52  receiver  will  be  discussed  first,  followed 
by  a  description  of  the  changes  caused  by  refueling  with  a 
KC-10  tanker. 

A  good  starting  point  for  the  air  refueling  process  is 
when  the  B-52  receiver  is  one  mile  behind  and  1000  feet 
below  the  tanker.  The  B-52  then  climbs  to  a  position  1/2 
mile  behind  and  500  feet  below  the  KC-135.  The  B-52  con¬ 
tinues  climbing  and  closing  on  the  tanker  to  a  point  slightly 
behind  and  below  the  tanker.  This  point  is  called  the  pre¬ 
contact  position  (approximately  fifty  feet  behind  and  ten 
feet  below  the  contact  or  refueling  position)  and  is  very 
important.  The  B-52  is  stabilized  (closure  rate  is  reduced 
to  zero)  at  the  precontact  position.  Precontact  is  impor¬ 
tant  for  two  reasons: 

1.  The  power  setting,  trim  (angle  of  attack),  and 
airspeed  are  approximately  the  same  for  the 
oontaot  position  (when  the  receiver  can  onload 
fuel ) . 

2.  The  aerodynamic  effects  between  aircraft  beoome 
much  more  notioeable  to  the  pilots. 
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As  the  B-52  moves  from  preoontaot  to  the  contact 
position*  the  reoeiver  pilot  attempts  to  keep  the  B-52 
on  tanker  center  line*  while  moving  the  aircraft  forward 
and  slightly  up.  This  is  often  a  difficult  task  for  the 
refueling  student  pilot  due  to  the  complexity  of  the  com¬ 
bined  effects  of  flying  two  large  aircraft  in  the  same 
airspace.  During  the  closure  to  the  contact  position* 
the  RC-135  must  put  in  nose-up  trim  to  maintain  level 
flight  conditions.  Approximately  seven  feet  from  the 
contact  position,  the  B-52  experiences  a  drag  increase 
(may  require  &  slight  throttle  increase),  and  a  nose  up 
trim  requirement.  However,  after  passing  that  point, 
there  is  a  slight  nose  down  trim  requirement  and  a  drag 
decrease  (may  require  reduction  of  throttle  setting). 

While  in  contact  with  the  tanker,  the  B-52  can  move 
laterally  four  degrees  (measured  from  the  tanker  boom) 
left  or  right  of  tanker  centerline  with  little  effect  on 
the  B-52  rolling  moment  (7:1-7).  If  the  B-52  moves 
laterally  more  than  four  degrees,  a  rolling  moment  is 
induced  on  the  B-52,  causing  it  to  roll  toward  tanker 
centerline. 

Refueling  with  the  KC-10  tanker  is  very  similar  to 
the  KC-135,  with  the  following  exceptions.  The  KC-10  has 
a  larger  wing  span  than  the  KC-135  (165  ft  vs  131  ft)  and 
also  <.  much  heavier  gross  weight  (590,000  lbs  vs  297,000  lbs). 
When  a  B-52  moves  from  preoontaot  to  contact  on  the  KC-10, 
there  is  no  increase  in  drag  at  the  seven  foot  point  as  with 
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th®  KC-135.  Also,  th®  B-52  oan  novo  ®ight  to  ten  degrees 
loft  or  right  of  tanker  centerline  without  getting  tho 
tondonoy  to  roll.  Finally *  tho  KC-10  boon  is  nuoh  longer 
then  that  of  the  KC-135  boon  and  keeps  the  receiver  farther 
behind  and  below  the  KC-tO. 

A  well  designed  sinulator  should  "fly"  the  sane  as  the 
aircraft.  Therefore*  the  previous  Info  mat  Ion  will  be 
qualitatively  used  to  oheok  the  results.  For  nore  Inf or- 
nation  on  B-52  refueling  procedures*  see  B-52  T.O. 

1 -1C-1  -1 5  (7). 
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III.  Analytical  Approach 


Mathematically  modeling  the  velocity  Induced  on  the 
reoelver  due  to  the  presence  of  a  tanker  was  done  by  three 
methods.  A  linear  circulation  distribution,  located  at 
the  quarter  chord,  vas  assumed  for  all  three  methods,  so 
that  the  equations  could  be  analytically  integrated  (See 
Appendix  A).  Methods  1  and  2  used  a  lifting  line  at  the 
quarter  chord  followed  by  vortex  sheet  to  represent  the 
tanker.  These  methods  did  not  simulate  roll-up  of  the 
vortex  sheet.  Method  3  used  one  horseshoe  vortex  with 
the  lifting  line  located  at  the  quarter  chord  and  two 
trailing  vortices  to  represent  the  tanker.  This  simulated 
two  fully  rolled  up  tip  vortices.  The  above  methods  were 
the  extreme  cases  (no  roll-up  and  fully  rolled  up  vortex 
sheets)  of  what  actually  occurs. 

The  linear  circulation  distributions  for  the  tanker 
in  Methods  1,  2,  and  3  were  selected  for  two  reasons. 

One,  to  keep  the  integration  as  simple  as  possible  and 
still  give  qualitatively  useful  results.  Two,  evaluate 
how  a  substantial  change  in  lift  distribution  affected 
the  prediction  of  induced  velocity. 

The  effeot  of  the  rolling  up  prooees  will  be  slight. 
Vortices  for  a  sweptwing  airoraft  are  fully  rolled  up  at 
three  to  four  span  widths  (11*13).  For  a  KC-10,  that 
amounts  to  660  feet  behind  the  wing  tip.  In  a  typical 
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situation,  a  raoaivar  is  only  25  faat  bahind  tha  tankar. 

Tha  dasoant  valooity  of  tha  vortex  pair  was  found  by  (11 *8) t 


vhara  b'  is  tha  distanoa  separating  tha  two  vortices 

(b'  *  trb/4  for  an  alliptioal  lift  distribution.)  For  a 
typical  refusing,  *  5.7  +  1.42  ft/sac. 

Notioa  in  Fig  4,  tha  B-52  flight  path  is  vail  balov 
tha  KC-10  tip  vortioea.  Evan  in  tha  contact  position, 
tha  raoaivar  is  twenty-five  faat  balov  tha  tankard  wing 
tips  and  tha  tip  vortices  are  only  beginning  to  roll-up, 
indicating  that  tha  affect  of  tha  rolling-up  process  should 
be  alniaal.  The  other  assuaptions  for  tha  analytical 
approach  are: 

1 .  Linearised  potential  flow  field 

2.  Raoaivar  is  on  tankar  oantarlina  (y«0) 

3.  Tankar  and  raoaivar  are  represented  by 
circulation  distributions  along  tha 
quarter  chords  (fuselage  effects  and 
tallplanes  ignored) 

4.  Swap  bn  ok  is  not  considered 

Both  -canker  and  raoaivar  vings  induced  a  down  wash  at  tha 
raoaivar  wings  quarter  chord.  Tha  total  dovnvash  (wj) 
on  tha  raoaivar  win*?  was  found  byt 

*2  *  V21  +  w22 

Wuera  w21  is  tha  lnduoad  valooity  oausad  by  tha  tankar  on 


the  receiver  end  v^z  ia  ^*e  velocity  the  receiver’ •  wing 
induces  on  itself. 

For  ell  methods  discussed  above,  equations  were  de¬ 
rived  to  express  the  induoed  velocity  at  any  point  in 
space.  To  get  the  tanker  influence  on  the  receiver 
noints  were  chosen  along  the  lifting  line  on  the  receiver's 
wing. 

1 

As  noted  above,  method  1  replaces  the  tanker  wing  by 
a  lifting  line  and  a  semi-infinite  vortex  sheet  combina¬ 
tion  (Fig  5).  Method  1  assumes  a  tent  shaped  circulation 
distribution  (Fig  1).  An  expression  for  the  velocity 
induoed  at  any  point  P  by  the  bound  vortex  (along  the  y 
axis)  and  the  planer  vortex  sheet  (in  the  x,  y  plane)  is 
derived  in  Appendix  A. 

The  velocity  induced  in  the  x  direction  (u),  is  much 
smaller  than  the  free  stream  velocity  (u«U)  so  it  is 
neglected.  The  velocity  induoed  in  the  y  direction  (v) 
will  always  average  to  be  aero,  provided  the  receiver  is 
on  tanker  centerline  which  was  an  assumption  for  the  study. 
Finally,  the  velocity  induoed  in  the  s  direction  (w),  is 
found  to  have  a  significant  affect  on  the  receiver's  flow 
field. 

The  induced  velooity  in  the  s  direction  (downwash)  is 
given  by  Eq  (22)  in  Appendix  A.  8q  (22)  was  nondimension- 
alised  by  dividing  by  the  free  stream  velooity  (U) .  Sub- 
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Vortex 


Sketch  of  Assumed  Tanker  Win; 


atituting  Eq  (32)  for  T+  yields  the  induoed  velocity  for 
any  point  in  spaoe  caused  by  a  tent  oiroulatlon  distribu¬ 
tion,  Thua: 


( 


21 


u 


4irb  L  (x2+i2)(R*b*/4+yb)*  (x2+s2)  (R+b2/4-yb)* 

'  -  kni-.-yV— 5  'i*  k.iiVjtfyrtoi 

b(x2+s2)R*  b  y2+a2+by+b2/4  b  a2+y2 

- 1*„  iliiai*  H| 

b  R+x  b  (R+ybtbZ/4)*+  x 

♦  1»B  3  (a) 

°  (R-yb-b2/4)*+  x 


In  Eq  (2)  the  following  information  la  needed  to  deter • 
mine  the  dovnvash  the  tanker  inducea  on  the  receiver: 

1.  Tanker  lift  ooeffioient  (C^) 

2.  Tanker  reference  area  (S) 

3.  Tanker  wing  apan  (b) 

4*  Point  of  intereat  P  (x,y,a) 

Method  2 

Similar  to  method  1 ,  method  2  alao  models  the  tanker 
wing  by  a  lifting  line  at  the  quarter  chord  and  a  semi- 
finite  vortex  sheet  combination  (Fig  5),  Method  2,  however, 
assumes  a  modified  linear  shaped  circulation  distribution 
(Fig  2),  The  only  difference  between  Method  1  and  Method  2 
is  the  shape  of  the  circulation  distribution)  therefore, 
the  discussion  for  Method  1  alao  applies  to  Method  2. 

From  Appendix  A,  the  induced  downvash  caused  by  a 
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modified  circulation  distribution  at  any  point  P(x»y,z) 
was  given  by  Eq  (25).  Equation  (25)  was.  also  nondimension- 
alized  by  dividing  by  the  free  stream  velocity.  By  sub¬ 
stituting  Eq  (33)  for  the  final  result  for  the  modified 
circulation  distribution  is: 


w 


21  = 


u 


“CLS 


-x(b/4+2y+4R/b)  .  x(bUbUR/b) 

+z2)(b2/l6+yb/2+R)4  (x*+z2)  (R+b2/4 


/4+yb) 


f 


,  .x(b-4yt4R/b?, .  .  . 

(x2+z2) (R-yb/2+b2/1 6)*  (x2+z2) (b2/4-yb+R)* 

2In t|tn|  r2+*2+b2M-yb  , 

b  y2+z2+b2/ 4+yb  b  y  2  +  z2+b2/ 1 6-yb/ 2 

h n |  tR.tgb^/J.<=  ).*■_- . * |  |  UlilW&A ) *  - 

b  (R+vb/2+b2/1 6)  ^  +  x  b  (Rs-v^b2/!)*  + 


-  X| 
X 


(R+yb/2+b  /  16)C  +  x  D  (RVy> '•b‘74)! 

+  4n|tR-?btb»/4)f  -  »|-gtn|(R-yb/2tb>/i6)t  - 

b  (R-vb+b2/4) *  +  x  b  (R-vb/2+b2/l6)*  + 


The  information  needed  to  determine  the  dc  mwash- the  tanker 
induces  on  the  receiver  in  Eq  (2)  also  applies  to  Eq  (3). 


Method  3 

The  previous  results  are  for  a  vortex  sheet  with  a 
linear  circulation  distribution,  simulating  no  vortex  roll 
up.  To  try  to  bracket  the  problem,  a  single  horseshoe 
vortex,  simulating  fully  rolled  up  tip  vortices,  may  be 
used  to  represent  the  tanker  wing.  The  derivation  for 
the  velocity  induced  by  any  horseshoe  vortex  is  found  in 
(1:191-201).  Figure  6-  shows  the  tanker  wing  represented 
by  a  horseshoe  vortex. 


Fig. 6.  Sketch  of  a  Horseshoe  Vortex  Modeling  the  Tanker 


Again,  select  point  P  on  the  lifting  line  of  the 
receiver  to  get  the  downwash  of  the  tanker  on  the  receiver. 
Then  vary  P  along  the  lifting  line  to  compute  the  downwash 
distribution. 

A  modified  te  wing  span  V  represents  the  effec¬ 
tive  wing  span  due  to  the  roll  up  of  the  tip  vortices. 

Muttray  (14)  states  the  correction  for  the  span  with 
a  constant  lift  distribution  is  .865.  therefore  b’^.SSSb 
was  used. 

From  the  3D  horseshoe  vortex  derivation  in  (1:197), 
the  induced  downwash  at  any  point  in  space  is: 

W21  = 


rH/  (x-x1 ) (y-y2)-(y«y1 ) (x-x2) 

(y-y-^ )  (z-22)-(2-z1  )  (y-y2))*+{  (x-x1 ) (z-z2)-(z-z1 )  (x-x2)}‘ 


+{ (x-x1 ) (y-y2)-(y-y1 ) (x-x2)} 

(x2-x1 )(x-x1 )+  (y2-y1 ) (y-y1)+(z2-z1 )(z-z1 ) 

{ (x-x.j  )2+(y-y1  )2+  (z-z1  )2}^ 


-(x2-x1 ) (x-x^J+Cyg-y^ ) (y-y2)+(82-a1 ) (z-z2) 
{ (x-x2)2+(y-y2)2+(z-z2)2}* 


♦. - ZCZ - - 

(z-z.,)  +  (y-y-| )' 


r  i  +  x-xi  _ l 

L  (x-x^ )2+(,j'-y-| )2+(z-z-j )2  J 


ty2-y) 

(z-z2)2+(y-y2)' 


[’  (*-*2)z+(y-y2)ii+(«-»2)i  ]  ) 


where  x,y,z  are  the  coordinates  of  point  P. 
As  indicated  in  Fig  6, 


Xl=x2=° 


zi=z2=0 


x  =  A 


z  =  H 


y1  =  -.865‘b/2 


,=  .86513/2 


Substituting  these  values  into  Sq  (4)  yields: 


=  ru/i-i^A  T.  865b(v+.  4325b)  __  -  .86?b(y-.4?2Sbj 


W01  =  1  H  I  I  .  I  Juft  I  .  y  ■  »  -  r  r.  ■■  ■■  ■  A  < 

*'  ^\bA2+H2b  l(.187b2+.'865yb+Rr  ( .1 87b<i-.865yb+R  ) 


+  (-.4325^7^ 

HZ+(y+.4325b)< 


[  (.isTeb^.sesyb+R)4] 


H2+(y-. 4325b)' 


P  +  - 
l  (.1 


87b2-.865yb+R)! 


The  circulation  is  found  from:  FH  =  ’  (s®e  Ap¬ 


pendix  B  for  derivation.) 


Dividing  Eq  (5)  by  U  yields  a  nondimens ional  equation 
for  downvash.  Thus,  Eq  (5)  becomes: 

V21  -  CLS  |  1.156AT  .865b(v+. 4325b) 

^  1  bA^+bH^  L  (.187b2+.865yb+R)* 

-  .865b(y-» 4325b)  1+  , l-.v412gb^) 

( ,187b2-.  86  5yb+R2)*-*  H2+(y+. 4325b)2 

1  +  A  1-  (,4125^) 

C.mb^+.SbSyb+RrJ  H2+(y-.  4325b)2 

1  +  _  A _ _  ]  I 

(.187b2-.865yb+R)* J  i  (6) 

In  Eqs  (2),  (3),  and  (6),  downvash  information  for 
the  receiver  flowfield  is  determined  by  varying  P  (x,y,z) 
along  the  receiver’s  lifting  line. 

Once  the  downvash  velocity  the  tanker  induces  on  the 
receiver  is  calculated,  the  relative  velocity  for  the 
receiver  can  be  determined.  This  relative  velocity  is 
found  by  the  vector  addition  of  the  free  stream  velocity 

and  the  tanker  induced  downvash  velocity.  The  aerodynamic 

* 

coefficients  are  calculated  using  both  the  relative  and 
free  stream  velocities,  and  the  differences  between  the 
coefficients  are  the  aerodynamic  interference  effects. 

In  the  above  equations,  changing  the  tanker  gross 
weight  is  accomplished  by  varying  the  tanker  coefficient 
of  lift.  Finally,  any  tanker  oan  be  modeled  by  changing 
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the  values  of  winespan  and  wing  area 


Compressibility  Considerations 

The  compressibility  effects  on  downwaeh  were  taken 
into  account  by  the  Goethert  transformation  (2).  The 
parameters  were  transferred  as  follows: 


Incompressible  SglBgtfUSll 


x 

y 

z 

b 


bcl 

x/B 

y 

z 

b 


where  B  -  (1+^)* 

A  discussion  of  the  results  obtained  using  Eqs  (2), 
(3),  and  (6)  and  comparison  with  other  data,  are  presented 
in  Chapter  Five. 


To  treat  more  complicated  configuration*  (sveepback, 
tailplana  af facts,  ato.)  a  nuaarioal  approach  is  required. 
The  vortex  lattice  method  (VLM)  vas  the  numerical  technique 
selected  for  study  because  of  its  well  established  and 
proven  ability  to  yield  aoourate  results  for  interfering 
lifting  surfaces  (15:89). 

In  Chapter  three,  a  single  horseshoe  vortex  vas  used 
to  represent  the  tanker's  wing.  In  contrast,  the  VLM  uses 
multiple  horseshoe  type  vortices  to  represent  a  wing.  The 
wing  is  overlaid  by  a  geometric  network  of  lattice  boxes 
or  panels.  A  typical  panel  contains  a  bound  vortex  and 
two  trailing  vortices.  The  bound  vortex  is  located  on 
the  panel's  quarter  chord.  A  panel  control  point  is  at 
the  center  of  the  three  quarter  chord  of  the  panel.  It 
is  at  the  control  point  that  the  flov  tangency  boundary 
condition  is  enforced.  See  Fig  7  for  a  typically  panelled 
win,?. 


Tig.  7,  Vortex  Lattioe  Panelled  Wing 


Pour  panels  vara  used  to  represent  tha  above  wing. 

Mora  aoouraoy  and  ohordvisa  raaolution  of  aarodynaalo 
foroas  oan  ba  aohiavad  by  using  ohordvisa  as  vail  as 
spanviae  panala.  To  aeoount  for  tip  affaotst  tha  last 
panal  vas  slightly  Indantad  fron  tha  ving  tip.  Tha 
aquation  usad  to  datarnlna  this  dlstanoa  Is  (15)* 

T«»  ■  (7) 

vhare  Mg  la  tha  nunbar  of  sanlspan  panels:  In  tha  abcva 
example,  Ng  la  4*  Tha  alight  lndantatlon  allova  tha  uaa 
of  laaa  panala  to  aehlava  tha  aama  accuracy  (15). 

Each  vortax  panal  lnduoaa  a  velocity  on  itaalf  and 
on  all  the  other  panala.  This  la  rapraaantad  by  an 
Influence  coefficient  matrix  A  where s 

tA,rj  *  (7.) 

A  typical  element  of  {A},  Ai j  represents  the  downwash 
velocity  at  the  control  point  of  panal  1  induced  by  a 
horseshoe  vortax  at  panel  j.  Fj  Is  the  circulation  at 
panel  j  and  la  the  flow  tangency  boundary  condition 
at  panel  j.  For  a  planer  ving  with  no  dihedral,  the  flow 
tangency  condition  is  "o  ■  -ua.  Gamma  (tha  circulation 
vector)  ia  than  found  by* 

F  »  (A}'1  o  (8) 
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Onoe  Otau  is  known,  the  aerodynamio  foree  distribution  eaa 
be  determined  (Saa  Appendix  D).  Rafaranoa  (1)  provides 
aora  datailad  daaoription  of  tha  vortax  lattice  method. 


( 


i 


The  following  asauaptiona  wara  aada  in  tha  vortex 
lattioa  approach! 

1.  Linaar  potential  flow 

2.  Winga  ara  rapraaantad  by  flat  plato a  (no  oaab 
or  twist) 

3.  Fuselage  affaota  ara  ignorad 

4*  Wing  tip  vortax  roll  up  ia  Ignorad 

5.  Raoaivar  ia  on  tankar  oantarlina  (y^  »0) 

I* 


An  existing  computer  program  called  LATEX  (Appandix  C) 
waa  aalactad  for  evaluation.  LATEX  uaaa  tha  vortax  lattioa 
method.  The  aignificant  advantagaa  of  LATEX  ara  ita  aim* 
plioity,  ability  to  handle  aora  than  four  lifting  aurfaoaa, 
and  direct  printout  of  induced  normal  velocity. 

The  coordinate  ayatam  waa  aat  up  auch  that  tha  a 
direction  ia  perpendicular  to  tha  tanker  wing  at  saro 
angle  of  attaok  (poaitiva  up),  x  la  in  tha  downatraam 
direction,  and  y  makes  an  orthogonal  system  going  out 
in  tha  direction  of  tha  tankers  right  wing.  The  coordinate 
ayatam  la  independent  of  tha  angle  of  attack,  therefore, 
to  use  tha  coefficient  in  tha  normal  direction  (Cs)  aa  the 
lift  coefficient,  tha  amall  angle  assumption  ia  made.  Tha 
origin  used  ia  always  at  tha  tankar  wing's  leading  edge 
centerline  (Saa  Fig  8). 

Two  methods  may  be  used  to  change  angle  of  attaok  in 
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the  LATEX  Program.  A  geometrlo  angle  of  attaok  cnerated  a 

| _ !  constant  o  (in  degrees)  for  the  antira  ayataa  of  lifting 

aurfaeaa.  A  aaoond  aathod  apaoiflaa  angle  of  attaok  (in 
radians),  by  aaana  of  a  table,  Tha  tabla  ia  part  of  tha 
input  fila  and  ia  uaad  to  ehanga  a  control  point  borundary 
condition  (See  Appendix  C).  Tha  advantage  of  using  tha 
input  or  tviat  tabla  ia  that  aaoh  wing  oan  now  ba  givan  a 
diffarant  angla  of  attaok.  Compariaona  uaing  both  tha 
gaoaatrio  a  and  tha  input  tabla  wara  aada  and  tha  rasulta 
indicated,  aa  axpectad,  that  aithar  aathod  provldaa  tha 
aaaa  results  for  tha  aaall  angles  considered. 

The  wake  was  modeled  to  extend  in  the  x  (downstream) 
direction.  In  the  caaa  of  the  KC-10  wing,  the  wake  extends 
at  i  *  0  at  a  constant  rata  equal  to  U.  Tha  wake  does  not 
roll  up  or  descend  as  an  actual  wake  would  do. 

Tha  majority  of  tasting  was  dona  with  four  lifting 
surfaces.  These  are  tha  KC-10  wing  and  tailplane,  and 
tha  B-52  wing  and  tailplane  (Saa  Pig  8).  Ten  aamispan 
panels  wara  uaad  for  tha  KC-10  wing,  twenty  semispan  panels 
were  uaad  for  the  B-52  wing,  and  fire  aamispan  panels  ware 
used  for  both  Sailplanes.  LATEX  reflects  an  image  of  eaoh 
wing  for  calculation  purposes,  ao  input  of  only  the  semi- 
span  data  is  necessary. 

Presently,  LATEX  has  no  diraot  method  to  determine 
off  body  information.  It  was  discovered  that  panels  with 
large  aspect  ratio  (0.5  ft  width  by  10  ft  length)  induced 
\  an  insignificant  amount  of  downvash  on  themselves  regard- 
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leas  of  panel  angle  of  attaok.  Whan  thaaa  panala  vara  uaad 
for  off  body  information#  the  lnduoad  dovnvaah  at  thaaa 
panala  vaa  due  to  tha  lnduoad  velooity  of  tha  flow  fiald. 
Tharafora,  thla  typa  of  panal  vaa  uaad  for  all  off  body 
data . 

Tha  pltohlng  aoaant  opaffloiant  (C^)  vaa  dataralnad 
by  aatabliahing  tha  B-52  at  a  daelrad  angla  of  attaok  and 
Itarativaly  aalaoting  a  tallplana  angla  of  attaok  to  aaka 
CM  *  0.  Tha  wing  angla  of  attaok  aalaotad  vaa  a  ■  6.0°, 
which  put  the  B-52  at  a  typloal  refueling  groaa  weight. 

Tha  B-52  oantar  of  gravity  waa  aalaetad  to  ba  at  tha 
quarter  naan  aerodynaaie  ohord  baeauaa  aoat  refueling  is 
accomplished  with  the  CQ  near  tha  quarter  ohord  position. 
Tha  B-52  was  then  put  behind  tha  KC-10  at  multiple  loca¬ 
tions.  The  B-52  lift  vaa  kept  constant  by  iteratively 
changing  its  angle  of  attack  at  each  position  behind  the 
tanker.  waa  determined  by  comparing  the  CM  for  the 

B-52  alone  with  tha  for  tha  B-52  behind  the  tanker. 

Tha  results  of  tha  Vortex  Lattioe  analysis  are  found  in 
Chapter  Five,  and  a  daaoriptlon  on  the  use  of  LATEX  is 
found  in  Appendix  C. 


) 
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The  results  of  ths  analytical  approach  and  the  LATEX 
ooaputar  prograa  are  oompsred  vlth  Douglas  KC-10  data  found 
in  (5)*  Thasa  data  vere  generated  by  a  Welaainger  Coaputer 
Method. 

The  Veissinger  aethod  replaoes  the  wing  by  a  plate 
of  sero  thickness  that  has  the  planfora  and  twist  identical 
to  the  wing.  The  ohordwise  loading  distribution  is  assuaed 
to  be  concentrated  into  a  lifting  line  located  at  the 
quarter-chord  of  the  wing.  Chordwise  control  points  are 
chosen  to  be  at  the  three-quarter-chord  line  and  the  flow 
tangenoy  boundary  condition  is  applied  at  these  points 
(no  flow  is  allowed  through  the  plate  at  these  points). 

The  lifting  line  and  its  trailing  vortex  sheet  are  treated 
to  be  continuous.  Substituting  the  vortex  pattern  for  the 
wing  and  applying  the  boundary  conditions  allows  the  for¬ 
mation  of  a  set  of  simultaneous  equations*  which  can  then 
be  solved.  For  a  more  detailed  analysis  of  the  Weissinger 
method,  see  NACA  Report  No.  921  by  DeYoung  and  Harper  (3). 

Tanker  lift  distribution,  tailplane  effects,  fuselage 
effects,  and  depression  of  the  vortex  sheet  were  considered 
by  Douglas  in  generating  the  downwash  data.  No  wind  tunnel 
data  or  flight  test  data  were  available  for  comparison. 

The  analytical  results  will  first  be  disoussed  followed  by 
the  vortex  lattice  results. 
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WyUs*;  fiaaa&ft 

Equations  (2),  (3),  and  (6)  wara  programmed  into  an 
HP-35CV  hand  oaloulator.  Tha  oaloulator  took  approximately 
6  saeonds  to  oompute  a  valua  of  downwash  at  aaoh  point. 
Recall,  that  tha  analytioal  solutions  raprasant  tha  antira 
tanker  by  a  rectangular  wing  with  a  linear  oiroulation 
distribution.  Downwash  vs  span-wise  location  on  tha 
receiver  wing  is  plotted  on  Figs  9-12. 

Figures  9-11  are  for  the  contact  position  (25  ft  below. 
125  ft  behind  the  tanker  leading  edge.)  Notice  on  Fig  9. 
the  tent  distribution  over-predicts  (in  comparison  with 
the  Douglas  data)  the  Induced  downwash  by  55$  at  a  span 
position  of  0.  The  modified  linear  and  the  horseshoe 
vortex  models  (with  b*  =  .865  in  the  latter  model)  agree 
on  the  average  within  25-35)1  of  the  Dbuglas  data.  In 
Fig  12  the  downwash  is  shown  for  the  precontact  position 
(35  ft  below,  175  ft  behind  the  tanker).  The  analytical 
curves  are  closer  to  the  Douglas  curve  at  a  distance 
farther  beind  and  below  the  tanker,  indicating  that  a  more 
accurate  model  is  needed  for  close  proximity  effects. 
Including  fuselage,  sweepback  and  tanker  tailplane 
influence. 

From  Figs  9-11 ,  it  can- be  seen  that  the  analytical 
equations  overestimate  the  average  downwash  over  a  90  ft 
semi-span.  The  modified  linear  distribution  gives  a  25$ 
disagreement  followed  by  the  horseshoe  vortex  with  a  35$ 
disagreement.  It  is  evident  from  these  figures  that  the 
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Induced  Downwash(??) 


KC-10  Data 

Weight  260,000 
Me0.6  - 

Sref*3647Ft 

Douglas  Prediction 

Tent 

Modified 

O - O 

Horse  shoe 

□ - □ 


Spanwise  Variation (Ft) 


Pig. 10.  Analytical  Downwash  vs  Span  Location 
(Contact  Position)  Wts26o,000  lbs 


Induced  Downwash 


S panwise  Variation (Ft) 


Fig. 11.  Analytical  Downwash  vs  Span  Location 
(Contact  Position) 
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Pig. 12.  Analytical  Downwash  vs  Span  Location 
(Pracontact  Position) 


shape  of  the  tanker  lift  distribution  and  sveepbaok  have 
a  substantial  effect  on  the  induced  downwash. 

While  the  analytical  equations  give  a  qualitative 
idea  of  the  tanker  flowfield,  better  results  can  be 
achieved  using  a  more  sophisticated  model.  The  vortex 
lattice  method  not  only  yields  better  results,  it  can  also 

evaluate  all  of  the  desired  interference  effects* 

\ 

VLM  Results 

The  vortex  lattice  method  was  used  to  study  the 
effects  of  wing  sweep,  tailplane  influence,  change  of 
bomber  pitching  moment,  and  how  the  presence  of  one 
aircraft  affects  the  lift  of  the  other  aircraft. 

The  vs  a  curves  are  shown  for  both  the  KC-10  and 
the  B-52  in  Figs  13  and  14*  Both  flight  test  derived  CL 
vs  a  curves  ere  greater  than  the  curves  generated  by  the 
vortex  lattice  method.  This  is  to  be  expected  since  the 
aircraft  wings  are  represented  by  flat  plates  in  LATEX 
(i.e.,  no  camber  or  twist)  where  in  reality,  both  wings 
have  camber  and  twist. 

Fig  1 5  shows  the  effect  on  induced  downvash  caused 
by  using  a  rectangular  wing  (with  the  same  span  and  wing 
area)  and  the  actual  swept  wing.  Results  are  plotted  for 
the  contact  and  the  precontact  position.  Notice  how  the 
swept  wing  induces  more  downwash  closer  to  centerline  and 
less  downwash  farther  away  from  the  centerline  than  the 
rectangular  wing.  The  figure  indicates  there  is  a  sig- 
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Pig. 13.  B-52  CL  vs  ot 


Pig.  14.  KC-10  CL  vs  o< 
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Fl«.  15.  VIM  Downwuh  Comparison  (at  ths  B-52  wing)  of 
Rsctangular  and  Swspt  Tanksr  Wines 


nifioant  difference  between  tha  swept  and  rectangular 
rasultai  therefore,  for  all  ooaputar  ganaratod  dovnvash, 
tha  actual  avapt  ving  planforaa  will  ba  uaad  for  greater 
aoouraoy. 

Sactlonal  lift  coefficient  (C, )  va  apan  for  tha  KC-10 
ving  ia  found  on  Fig  16.  Tha  ving  vaa  panelled  vlth  one 
chordviaa  panel  and  tan  aeaiapan  panala.  Fifty  panala  and 
100  panala  par  aeaiapan  vara  alao  evaluated.  It  vaa  found 
that  vhan  one  chordviaa  panel  and  tan  aeaiapan  panala 
represented  the  tanker  ving.  an  increaaed  nuaber  of  panels 
gave  no  significant  improvement  (less  than  Iff  change  in 
induced  dovnvash). 

The  dovnvash  field  behind  the  KC-10  was  evaluated 
vith  and  without  the  tanker  tailplane.  Figures  17-20 
indicate  the  induced  dovnvash  vs  span  froa  h  alia  to  the 
contact  position. 

Notice  there  is  only  a  slight  Increase  in  induced 
dovnvash.  caused  by  the  tanker,  as  tha  receiver  aovea 
up  toward  the  precontact  position.  At  the  contact  and 
precontact  positions,  the  results  are  compared  vith  the 
Douglas  data.  Figures  19  and  20  indicate  that  tha 
induced  dovnvash  is  slightly  in  error  vith  or  vithout 
tha  tailplane  vhan  tha  KC-10  angle  of  attack  ia  6.8 
degrees.  HovaVtr,  Tig  21  indicates  that  tha  KC-10  vith 
tha  *  ilplane  aatchaa  tha  Douglas  data  for  an  angle  of 
aliack  of  2.3  dagraaa.  Therefore,  inclusion  of  tha 
tanker  tailplane  results  in  aora  accurate  values  of 
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induced  downwash,  especially  at  a  low  KC-10  angle  of  attack. 
Also,  tailplane  effects  are  far  less  critical  at  points 
farther  behind  and  below  the  KC-10. 

The  tailplane  must  be . included  if  tanker  pitching 
moments  are  required.  Figure  22  indicates  that  the  B-52 
induces  a  significantly  larger  upwa^h  on  the  KC-10  tail- 
plane  than  on  the  wing.  This  causes  a  nose-down  pitching 
moment  on  the  tanker  supporting  the  need  for  the  tanker 
to  trim  nose  up  when  the  bomber  moves  from  precontact 
to  contact  position. 

When  the  B-52  flies  behind  the  KC-10,  the  change 
in  upwa3h  in  the  tanker  flowfield  also  causes  a  change  in 
the  KC-10  lift  coefficient.  This  change  is  shown  in 
Fig  23  at  various  distances  behind  and  below  the  tanker. 
There  is  a  definite  increase  in  lift  coefficient  at  all 
vertical  distances  when  the  B-52  is  125  ft  behind  the 
KC-10' a  wing  leading  edge.  However,  the  maximum  change 
of  tanker  C^  is  only  2.92$  of  the  value  without  the 
B-52,  indicating  that  the  bomber  wing  has  a  minimal 
effect  on  the  tanker’s  wing. 

Considering  the  lift  increase  on  the  tanker  wing  by 
itself,  the  pilot’s  reaction  is  to.  slightly  decrease  the 
angle  of  attack  to  maintain  level  flight.  However, 
coupled  with  the  reaction  of  the  tanker  tailplane  due  to 
the  receiver,  the  overall  effect  is. a  need  for  noseup 
trim  to  maintain  level  flight.  The  final  tanker  angle 
of  attack  is  slightly  less  than  the  original  (due  to  the 
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lift  increase)  and  the  need  for  hoseup  trim  is  due  to  the 
downward  pitching  caused  by  the  receiver. 

The  LATEX  Program  was  run  for  numerous  positions 
behind  and  below  the  tanker..  Figure  24  indicates  the 
points  checked  and  the  dark  black  line  is  the  B-52’s 
desired  flight  path.  The  object  was  to  determine  the 
B-52’s  change  in  angle  of  attack  when  approaching  the 
KC-10  and,  directly  related  to  the  angle  of  attack,  the 
change  in  the  B-52  pitching  moment.  The  results  are 
shown  on  Figs  25  and  26. 

Both  pitching  moment  and  angle  of  attack  are  sub¬ 
stantially  effected  by  the  distance  below  the  tanker 
(approximately  50%  decrease  for  a  50  ft  decrease  in 
receiver  altitude).  This  result  was  al3o  expected  and  is 
justified  by  both  the  Douglas  data  and  flight  experience. 
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VI.  Conclusions  and  Recommendations 


The  idee  that  a  tanker  can  be  represented  by  a  linear 
circulation  distribution  and  a  vortex  sheet  or  a  single 
horseshoe  vortex  proved  to  be  reasonable.  Equations  (2), 
(3),  and  (6)  provide  qualitatively  accurate  results  for 
induced  dovnvash  velocity  at  the  reoeiver  and  can  be 
easily  programmed  into  a  hand  held  programmable  calculator. 
This  method  is  sufficient  for  only  a  few  data  points  and 
is  a  good  place  to  start  to  get  an  approximate  idea  of  the 
flow  field.  It  also  shows  that  the  changes  in  the  tanker's 
lift  distribution  significantly  influences  the  induced 
downwash. .  The  analytical  solution  generally  over-estimates 
the  induced  downwash  by  25— 35% *  because  the  analytic  solu¬ 
tion  uses  a  rectangular  wing,  and  the  downwash  fields  of 
rectangular  and  swept  wings  differ  (8). 

The  Vortex  Lattice  approach  provides  a  more  accurate 
method  that  is  capable  of  predicting  not  only  the  induced 
downwash,  but  also  the  change  in  receiver  angle  of  attack, 
change  of  the  receiver  pitching  moment,  and  change  in  the 
tanker  lift  coefficient. 

Other  areas  of  interest  were  also  studied.  The  com¬ 
parison  of  rectangular  vs  swept  wings  indicated  that  the 
actual  swept  wing  model  provides  a  closer  estimation  of 
induced  downwash  and  should  be  used  in  air  refueling 
evaluation.  The  overall  effect  of  the  tanker  tailplane 


49 


on  the  bomber  was  found  to  be  relatively  small  (less  than 
10$  change  in  Induced  dovnwash),  especially  with  the 
receiver  veil  below  the  tanker. 

Comparison  of  the  results  of  LATEX  and  Douglas  data 
was  found  to  be  very  good  (See  Figs  19*  20,  21),  especially 
at  lover  tanker  angles  of  attack.  The  trends  and  magni¬ 
tudes  of  the  change  in  the  receiver's  angle  of  attack  and 
pitching  moment,  and  change  in  tanker  lift,  were  in  the 
proper  direction  and  reasonable  based  on  flight  experience. 

LATEX  Is  capable  of  matching  Douglas  predictions  of 
induced  dovnwash  within  5-10$,  by  representing  the  tanker 
and  receiver  aircraft  by  flat  plate  planforms.  Results 
indicate  that  the  VLM  may  provide  the  majority  of  infor¬ 
mation  necessary  for  evaluating  air  refueling  aerodynamics. 
However,  the  total  error  caused  by  using  VLM  or  any  computer 
prediction  will  not  be  realized  without  wind  tunnel  and/or 
flight  test  data, 

A  number  of  assumptions  were  made  in  this  study. 

These  include  neglecting  effects  of  the  fuselages,  wing 
camber,  and  wake  roll  up.  Effects  of  out  of  symmetry 
plane  approaches  on  the  receiver  rolling  moment  were  also 
not  treated. 

Fuselage  effects  may  have  a  significant  influence 
on  the  pltohing  moments  of  both  the  receiver  and  the 
tanker.  The  fuselage  may  also  vary  the  dovnwash  flovfleld 
by  interrupting  the  dovnwash  caused  by  the  wings.  The 
results  of  a  computer  analysis  with  and  without  the  tanker 


fuselage  should  be  compared  with  the  results  obtained  here 
to  determine  the  importance  of  modeling  the  fuselage. 

Wind  tunnel  tests  would  also  prove  beneficial.  They  would 
give  more  credibility  to  both  the  contractor  data  and  the 
LATEX  results  and  also  give  a  better  idea  of  what  aotually 
happens  to  the  flow  behind  the  tanker. 

Another  area  of  investigation  is  to  study  the  effeots 
of  camber.  The  LATEX  program  needs  to  be  modified  to 
permit  treatment  of  cambered  wings.  This  would  determine 
if  a  flat  plate  wing  is  a  good  assumption. 

A  more  accurate  representation  of  the  tanker's  wake 
could  also  be  modeled.  Wind  tunnel  tests  are  necessary 
to  determine  the  location  of  the  wake,  then  a  computer 
model  could  be  developed  to  simulate  the  wake. 

A  final  area  to  investigate  is  the  receiver  rolling 
moment  changes  when  approaching  the  tanker.  This  area  is 
critical  for  a  realistic  simulation,  because  it  is  a  prob¬ 
lem  area  for  many  student  pilots.  It  may  be  possible  to 
modify  LATEX  to  provide  the  information  that  is  needed. 

By  accomplishing  these  studies,  the  total  potential 
of  using  a  vortex  lattice  method  to  model  air  refueling 
between  a  tanker  and  heavy  receiver  would  be  realized. 
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)  Appendix  A 

Velocity  Field  Induced  by  a  Linear  Lifting  Line 
and  a  Semi-Infinite  Sheet  at  Any  Point  in  Space 


bound 

vortex 


z 


Fig.  27.  Sketch  of  Assumed  Tanker 
The  velocity  induced  at  point  P  is  caused  by  the 
bound  vortex  (lifting  line)  located  on  the  y  axis  and 
a  planar  vortex  sheet  issuing  from  the  bound  vortex. 

The  sheet  is  located  in  the  x-y  plane  and  all  vortex 
legs  are  parallel  to  the  x  axis.  The  effect  due  to  the 
bound  vortex  will  be  considered  first. 

The  induced  velocity  (V)  at  P  due  to  a  vortex  fila¬ 
ment  segment  dl  from  (1)  is: 

dVu  =  rnj?1 


(9) 


where  r  is  a  vector  from  dl  to  P.  The  ££  subscript  refers 
to  the  lifting  line. 

In  this  case: 

d7  ■  dyj  r  =  xT  +  (y-y^)J  + 


therefore; 


d  £  x  dr  *  -Txdy1  +  Tzdy1 

r3  =  {x2+(y-y1)2+z2}3^2 


x\y\z^  indicates  the  coordinate  on  the  vortex 
Substituting  in  to  (9): 

dV**  *  ^  r££^yi|^ady1  -  *  ^(y1)^1 

4tt{  x2+  ( y-y^+z2}  4ir{x^+(y-y  1)2+z2} 3/2 


(9a) 


To  find  the  total  contribution  of  the  lifting  line, 

Eq  (9a)  must  be  integrated  over  the  wing  span,  therefore: 
b/2  1,1  b/2  .  i 

dvcff  *  T  t  r££^y  ^zdy  -  I  f  _r££  7  ^Xdy 

{x^(y-y1  )2+zV'2  ^V(x2+(y-y1  )2+z2}3^2  (10) 

-b/2  -b/2 

Before  proceeding,  the  contribution  of  the  vortex 

sheet  will  be  derived  in  integral  form. 

The  vortex  sheet  consists  of  an  infinite  number  of 

filaments  eminating  from  the  lifting  line  and  extending  to 

infinity.  The  strength  of  each  filament  is  ^..(y1) 

—  dy ' 
dy1 
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in  accordance  with  Helmholtz 1  Vortex  theorem  (6:148) » 

For  one  filament  issuing  from  y  on  the  lifting  line* 
the  induced  velocity  at  P  due  to  this  filament  is 

av,»  -drtt(y1)iy1  hxt  (id 

dy  4rr3 

dT  is  now  a  small  segment  on  the  filament  and  again  r  is 
a  vector  from  dT  to  P. 

In  this  case: 

r  =  (x-x^T  +  (y-y^T  + 
dT  *  dx1! 


Therefore: 


dTxr  *  Tdx1 (y-y1 J-Jzdx1 

r3  =  { (x-x1 )2+(y-y1 )2+a2)3^2 
Substituting  into  Eq  (11)  gives: 

av0  .  d.ru.(^2^1  .  *W.-W  (v-yl) 

dy1  4ir  ( (x-it1  J^Ky-y1  )!+zZ)3^2 

The  contribution  of  the  entire  filament  is  found  by  in¬ 
tegrating  Eq  (12)  over  x1  giving 


(12) 


"■  '  8' 


tTa-Kty-y1 )  dxl 


d7  ^  ‘  {"(x-xf)  +  (y-y1)  +  «2>3/2 


(13) 


x1  »0 


Integrating  Eq  (13)  over  the  span  gives  the  equation 

Induced  velocity  due  to  the  vortex  sheet.  Therefore, 
b/2  «o  . 

?  =  i  )  +T*-K(y-y  )  i  i 

'  *7/^7 ((x-xV^y-yW.3'*  dXdy  (H) 


-b/2  o 


The  total  velocity  contribution  is  found  by  adding  the 
results  of  the  sheet  and  the  lifting  together  (Eq  (10) 


plus  Eq  (14))  giving: 

b/2 


(X.y.*)  iltj  {x*+(y-y 
-b/2 
b/2  » 


♦  k  f  f~~\ - dxV 

**  J  /dy1((x-x1);!+(y-y1)  +*  >3/ 


-b/2  o 
b/2  . 


'hlf  f  - dr^5)l'if'yi)5  5  jy;  dxV 

**  J  J  dy1{(x-x1)2+(y-y1 )2+a2}372 


-b/2  o 


+  f—s — -(4v"  2  in  dy1f  (15> 

J  { x2+ (y-y^ )2+z2} 3' 2 
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Intagrating  over  x1  yields  equations  of  the  following 
forint 


Integrating  and  'evaluating  at  the  limits  gives: 


(y-y1  r+ a  (x^+ty-y1  r+z2} *{  (y-y1  )2+z2} 


Substituting  back  into  Eq  (14)  yields: 
b/2 
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When  these  values  are  substituted  into  Eq  (16),  all 
integrals  can  be  put  in  standard  form  and  be  integrated. 

Substituting  the  linear  form  of  Gamma  into  Eq  (16) 
the  induced  velocity  in  the  x  direction  (u)  becomes 


.1/bv*  1 


Cl  -y1  /|)  dy1 


u  =  io*  r  v  r  v  y  *  (1?) 

^ J  {R+(y1)^-2yy1}3/2  J  { R+  (y1  J^yy1 } 3/2 


-b/2 


Integrating  Eq  17  yields: 


■b/2-2v+2R/b  +b/2+2y+2R/b  - 


4tt(x2+z2)  (R-yb+b2/4)4  (R+yb+b2A)’ 


} 

b(R)4 


(18) 


Again,  substituting  the  linear  form  of  Gamma  into 
Eq  (16)  the  induced  velocity  in  the  y  direction  (V  ) 
becomes: 


v  =  z  r  2T  o  dy 


2ro  dy 


^  J  b{ (y-y)^+z^}  J  b{ (y-y^ )^+z^} 


-b/2 


7 


2r  x  dv; 
o 


.f  _  2r0  X  dy1 

J b{x^+(y-y1 )2+z*> 


b{ x2 ' (y-y ' )2+z2} *{ (y-y ' )2+z2} 


-b/2 


!{(yV)V) 

(19) 


Integrating  Eq  (19)  yields: 
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V 


-  [  2  tan-”'  (y/z)  -  tan”"*  b/2-y  +  tan”3  b/2+4 

27Tb  z  z 


-  2  tan 


-1 


(x2+y^+z2)* 


-  tan' 


- — ,  . 

{ (y+b/2)z+x2+a2)® 


+  tan 


-1 


('(r-Va!^!4)*3 


(20) 


Finally,  to  find  the  induced  velocity  in  the  z 
direction  (downwaeh)  substituting  Gamma  into  Eq  (16)  gives: 


b/2 


w  =  - 


rl  rfl°  f  -  f£o  /*ixV,)dy 

4*  L^V  (yV)O  ^"7  (y-y1  )2+z2 


-b/2 


0+ 


+  fV/ - r-^.M  1  a  g 

b  J  (x2+(y-y1)2+z2}^{(y-y1)2+z2} 


-b/2 

b/2 


0- 


.  2v  r  (y-y1)^1  t  xr 
jTsry  ° 


7 


(y-y  )  +z‘ 


/ 


(l+yV^Jdy1 

{x2+(y-y3)2+z2}3/2 


-b/2 


b/2 


+  xr  /■  (i-yV^dy1  (2d 


( 
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-m 


“*3^ins?pjr" 


Integrating  Eq  (21  )j 


w 


2  2 
Z4+  v 

"2  ? 

I  1 _ Ll.*-  , 


z‘+y*+by+b*/4 


-  2 

b  (R)*  +x 


1 

¥ 


ln|iRiWZi4^| 

(R+yb+b*/4)*  +x 


+  1 
¥ 


In 


(R-yb+b2/4)*  -x  +  x(b/2+2y+2R/b) 
(R-yb+b^/4)4  +x  (x2+z^)(R+b2/4+yb)* 


+  x 


(b/2-2y±2R/b) 

..2,  2\/niv2/i 


(x2+.s2)(R+b2/4-yb)6 


;U|/b)  ]  (22) 

(xZ+*2)(R)* 


Equation  (22)  gives  the  induced  velocity  in  the  z  direction 
at  any  point  in  space  (except  x^y^z^O)  due  to  a  linear  tent 
circulation  distribution  and  its  vortex  sheet. 

To  determine  the  effect  of  tanker  lift  distribution 
on  induced  velocity*  a  modified  linear  circulation  dis¬ 
tribution  was  also  evaluated. 


z 


Fig.  29  Tanker  Represented  by  Modified  Linear 
Circulation  Distribution 
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The  only  difference  between  the  modified  linear 
circulation  distribution  and  the  tent  distribution  is  the 


value  for  Gamma. 
In  this  case: 


2r 


M  (1+V5) 


-b  <  y 

2 


1 


<  -b 

x 


r(y')  =  rM 


-b  < 

I 


<  b 


z 


b  <  v1  <  b 

Z  2 


and 


2rM 

*75 


_b  <  <  -b 

5  Z 


»b  <  y^ 

z 


<  b 


Z 


-2rM 

T75 


b  <  y 

z 


1 


<  b 
2 


Substituting  these  values  for  Gamma  into  Eq  (16) 
and  integrating  yields: 


u  = 


/M  -  [  •  =^2-y-,4R/bi  ♦ 

(x2+z2)4*  2{x2+(y+b/4)2+22)*  {x2+(y+b/2)2+z2}* 


+  ■b^.tym/b  +  .  .btAR/b  .......  ] 

2{x2+(y-b/4)2+zM  {x^y-b/^+z5}* 


(23) 
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Appendix  B 


Comparison  of  Elliptic  Circulation  Distribution  and 
Linear  Circulation  Distributions 

In  order  to  compare  the  results  of  the  various  tanker 
circulation  distributions  the  area  under  the  circulation 
vs  span  curves  was  made  equal.  This  in  effect,  caused 
the  total  lift  to  remain  constant,  even  though  the  dis¬ 
tribution  was  allowed  to  vary.  In  the  equation  for 
elliptical  circulation: 

r(y’)  -  re[  1.(^2]  i  -j  '  y1  <  b  (26) 


r 


e 


JUu 

pUb 


(27) 


where  r  is  constant  for  each  flight  condition. 
For  the  tent  distribution 


-b  <  y1  <  0 

2 


0  <  v^  <  b 
y  2 


(28) 


By  integrating  Eqs  26  and  28  over  the  wing  span  and 
equating  the  results,  a  value  for  can  be  determined. 
Equation  (26)  will  be  integrated  first. 
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for  an  elliptical  distribution. 

Now  integrate  Eq  (26)  for  the  tent  distribution 
b/2  0  b/2 

/r(y1)ay1  7'*  (1-&)*1 

-b/2  -b/2  0 

Tt  (b/2) 


s 


(3C 


Equating  Eqs  (29)  and  (30)  ensures  the  total  cir¬ 
culation  and  lift  is  constant 


Solving  for 


(31) 


Using  Eqs  (32),  (33).  and  (34)  for  functional  forms 
of  circulation,  caused  the  total  lift  to  be  constant. 
This  allows  direct  comparison  of  the  results  of  the 
equations  in  Chapter  III. 
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Appendix  C  . 


Operation  of  the  LATEX  Program 


Background 

The  LATEX  program  originated  at  Rockvell  International 
and  was  called  Unified  Vortex  Lattice.  It  was  modified 
repeatedly  by  ASD/XRHI  at  Wright-Patterson  AFB  for  their 
particular  needs.  One  of  these  modifications ,  called 
LATEX  IV,  was  used  for  evaluation  in  this  report.  LATEX 
was  run  on  a  VAX/UNIX  computer  system  at  the  Air  Force 
Institute  of  Technology.  The  average  turn  around  time 
(from  Input  to  Output)  was  about  2  minutes.  The  re¬ 
mainder  of  this  appendix  consists  of  a  brief  user  guide 
that  describes  how  to  build  an  input  file  and  also  how  to 
interpret  the  output  file. 

Data  Input 

Input  data  for  LATEX  is  accomplished  by  a  properly 
formatted  file.  Formatting  information  is  found  on  Figs 
30-32  and  a  sample  input  file  is  found  on  Fig  33.  The 
following  is  a  brief  discussion  on  how  to  set  up  an  input 
file. 

The  first  entry  contains  alphanumeric  title  infor¬ 
mation.  Data  input  starts  with  line  number  1  in  the  12th 
column.  Line  numbers  are  left  justified  and  go  in  columns 
10,  11,  and  12.  All  other  data  is  right  justified -and 


Pisr.31.  Formating  Information(cont. 


Fi<*. 32.  Formating  Informat i on (cont. 
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Fis. 33.  LATEX  Sample  Input  File 


starts  in  columns  14,  25,  37,  49,  and  61,  For  example  in 
Fig  33,  line  1;  starting  in  column  14  is  the  mach  number 
0.6,  starting  in  column  25  is  O.C  meaning  symmetric 
geometry,  0,0  starting  in  column  37  means  a  new  A  matrix 
is  calculated,  2, Cl  starting  in  column  49  prints  out  the 
information  needed  for  the  air  refueling  problem,  and 
column  61  is  not  used.  Line  15,  columns  37-40  have  a 
201 ,  indicating  the  twist  table  starts  at  line  201 .  The 
twist  table  can  be  used  to  change  the  angle  of  attack  of 
one  or  all  of  the  lifting  surfaces  and  will,  be  discussed 
later.  In  line  20,  the  4  indicates  there  are  4  lifting 
surfaces.  The  50,  100,  150,  175  are  the  relative  location 
of  the  data  for  each  lifting  surface.  In  this  case  the 
KC-10  wing  starts  at  50,  the  B-52  wing  starts  at  100,  the 
B-52  tailplane  starts  at  150,  and  the  KC-10  tailplane 
starts  at  175.  Lines  50-70  are  information  about  the 
KC-10  wing.  In  line  50  the  1  in  column  14  means  the  KC-10 
wing  is  lifting  surface.  The  0.0  in  columns  25-27  means 
a  symmetrical  image  is  calculated.  The  1.0  in  columns 
37-39  means  there  is  one  chordwise  panel  and  the  10.  in 
columns  49-51  means  there  are  10  span  wise  panels.  The 
3.0  in  columns  61-63  means  the  KC-10  wing  is  broken  up 
into  3  sections  or  elements  (Fig  34). 

The  twist  table  starts  at  201  and  goes  to  236  (or 
higher  if  required) .  The  twist  table  can  change  any  span 
wise  panel  angle  of  attack  by  changing  its  boundary  condi¬ 
tion.  Inputs  must  be  in  radians  and  a  minus  input  corre- 
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82.68* 


Panellina:  Used  For  KC-10  Wing 


sponds  to  a  positive  angle  of  attack.  In  Fig  33  lines  201 
and  206  are  all  zeros,  indicating  that  the  tvist  table  was 
not  used  for  the  KC-10  wing,  therefore,  it’s  angle  of 
attack  is  6.0°.  Lines  211-231  are  -.0285  radians  (1.63°) 
indicate  the  B-52  wing  angle  of  attack  is  7.63°  (1.63°  + 
6.0°).  The  B-52’s  tailplane  A0A  is  3.9°  (6j0  -  ( . 0367)  (180) 
and  the  KC-10* s  tailplane  A0A  is  3.0°  (6.0  -  ( . 0524) (1 80) . 
The  advantage  of  using  the  twist  table  is  that  all  four 
lifting  surfaces  can  be  placed  at  different  angles  of 
attack. 

The  last  line  of  data  must  contain  a  minus  sign  in 
column  1  to  indicate  the  end  of  input  information.  In 
Fig  33  this  occurs  on  line  236.  The  line  following  the 
last  data  line  is  another  alphanumeric  title  card.  The 
next  line  must  have  a  minus  sign  in  column  1  and  repeat 
some  information  found  in  the  data.  In  Fig  33  line  1 , 
Mach=0.6  is  repeated.  The  final  line  must  have  a  1  in 
columns  5  and  10. 

The  file  can  now  be  copied  onto  the  IN  file  and  LATEX 
can  be  run  for  the  desired  input  data. 

Data  Output 

The  output  for  the  sample  input  file  is  found  on 
Figs  35-38.  When  the  LATEX  is  run,  OUTPUT  is  the  file 
name  for  the  desired  output  data.  On  Fig  35  the  input 
information  is  printed.  The  values  of  XQ,  IQ,  and  ZQ  are 
the  x,y,z  location  of  each  panel  control  point  XQV  is  the 
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x  value  of  the  bound  vortex.  Delta  Cp  is  the  change  of 
center  of  pressure  for  the  back  panel.  Delta  CDIIC  is 
the  change  in  the  panel’s  induced  drag. 

The  induced  normal  wash  at  the  vortex  is  w__.  This 
value  is  negative  for  dovnwash.  Under  the  section  "spanwise 
loads,"  T  and  Z  are  the  y,z,  coordinates  of  the  loads, 
cx,  cy,  and  cz  are  s  action  coefficients  in  x,y,z  direction. 
When  the  small  angle  approximation  is  made,  cz*c£»  Under 
the  section  "panel  parameters"  in  Figure  38,  Panels 
1,  2,  and  3,  represent  the  3  elements  intheKC-10  wing, 
likewise  panels  4*  5  and  6  represent  the  3  elements  of  the 
B-52  wing,  panel  7  is  the  bomber  tailplane  and  panel  8  is 
the  tanker  tailplane.  To  find  for  the  bomber  wing,  the 
values  of  czp  listed  by  panel’s  4.  5  and  6  are  added  to¬ 
gether  then  multiplied  by  2  (due  to  symmetry).  The  total 
parameters  are  for  the  entire  system.  For  example,  the 
value  for  cz  is  found  by  adding  all  the  values  in  the  csp 
column  and  then  multiply  by  a  factor  of  2, 

The  previous  discussion,  along  with  the  sample  run, 
should  provide  the  information  necessary  to  use  LATEX. 
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Appendix  D 


Relationship  of  Aerodynamic  Coefficients 
and  Circulation  (VLM) 

As  previously  mentioned,  the  VLM  used  a  number  of 
discrete  spanwise  and  chordwise  panels  to  simulate  a 
wing.  For  each  panel,  a  circulation  (r)  was  determined. 
Once  T  was  known,  the  aerodynamic  coefficients  (C^,C^, 

C^)  were  found  as  shown  below. 

Section  Lift  Coefficient 

From  (1:113),  it  was  shown  that  the  section  lift  was 
found  by: 


Ct  =  (35) 

.  5pu  c 


c 

where  l  =  J"  pUy(s)ds  (36) 

o 


c 

For  the  VLM,  f  y(s)ds  *  T,  where  T  was  the  sum  of  all 
o 

Ps  in  the  chordwise  strip  of  interest.  For  example,  if 
there  were  2  chordwise  panels,  r  *  + Tg*  The  section 

lift  could  then  be  reduced  to: 


2.  =  pUT  (37) 


Substituting  Eq  (37)  for  l  in  Eq  (35),  the  section 
lift  coefficient  could  now  be  determined  with  known  infor- 
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matlon 


Moment  Coefficient 

The  moment  caused  by  aerodynamic  forces  acting  on  a 
lifting  surface  was  dependent  on  the  reference  axis. 
Specific  x,y,z  coordinates  of  the  reference  axis  were 
entered  into  the  input  file.  Once  the  reference  axis 
was  selected,  the  moment  coefficient  could  then  be  deter¬ 
mined  by: 
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(42) 


The  moment  (M)  was  found  by  multiplying  the  aerody¬ 
namic  force  by  the  distanoe  avay  from  the  reference  point. 
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